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ABSTRACT 

We present C/iandra observations of 17 optically-selected, X-ray weak narrow- 
line Seyfert 1 (NLSl) galaxies. These objects were optically identified by Williams 
et al. (2002) in the Sloan Digital Sky Survey Early Data Release, but were not 
found in the ROSAT All-Sky Survey (RASS) despite having optical properties 
similar to RASS-detected NLSls. All objects in this sample were detected by 
Chandra and exhibit a range of 0.5-2 keV photon indices F = 1.1 — 3.4. One 
object was not detected in the soft band, but has a best-fit F = 0.25 over the 
full 0.5-8 keV range. These photon indices extend to values far below what 
are normally observed in NLSls. A composite X-ray spectrum of the hardest 
objects in this sample does not show any signs of absorption, although the data 
do not prohibit one or two of the objects from being highly absorbed. We also 
find a strong correlation between F and Likev; this may be due to differences 
in Lboi/-^Edd among the NLSls in this sample. Such variations are seemingly in 
conflict with the current paradigm that NLSls accrete near the Eddington limit. 
Most importantly, this sample shows that strong, ultrasoft X-ray emission is not 
a universal characteristic of NLSls; in fact, a substantial number may exhibit 
weak and/or low-F X-ray emission. 

Subject headings: galaxies: Seyfert — galaxies: active — quasars: general — X-rays: 
galaxies 



1. Introduction 



Osterbrock & Pogge (1985) initially defined narrow-line Seyfert 1 galaxies (NLSls) by 
their striking optical spectral characteristics: strong, narrow H/? emission (later formally 
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defined to be FWHM < 2000 km s"^ by Goodrich 1989), weak [O III] relative to H/3, and 
strong Fc II. These properties put NLSls at one extreme end of the so-called Boroson & 
Green (1992) "Eigenvector 1," thought to correspond to emission from lower-mass nuclear 
black holes coupled with near-Eddington accretion rates (Boroson 2002). X-ray observations 
have revealed a strong soft X-ray "excess" in NLSls (e.g., Lcighly 1999), further bolstering 
the high Lboi/-^^Edd — low-mass black hole hypothesis (Pounds, Done, & Osborne 1995; Wang, 
Brinkmann & Bergeron 1996). Indeed, Boiler, Brandt, & Fink (1996, hereafter BBF96) found 
a possible anticorrelation between X-ray spectral slope and H/3 line width, with NLSls gener- 
ally having softer X-ray spectra than other AGN. Ultrasoft X-ray selection has consequently 
proven to be an essential tool for the discovery of large numbers of NLSls (e.g. Grupe 2000; 
Grupe et al. 2004). 

The disadvantage of selecting NLSls solely upon their X-ray properties is that it can in- 
troduce into NLSl samples a strong bias toward those exhibiting an ultrasoft excess (Forster 
1999). Since NLSls are primarily defined by their optical properties, the true nature of their 
X-ray emission is thus difficult to determine. Though it is well known that some NLSls are 
ultrasoft X-ray sources, previous samples of optically selected NLSls were simply too sparse 
to determine how many, as well as whether or not a significant number of NLSls have hard 
X-ray spectra. With the advent of the Sloan Digital Sky Survey (SDSS; York et al. 2000), 
it is possible to build large catalogues of NLSls with homogeneous selection criteria based 
on their optical spectra alone (see WiUiams, Pogge & Mathur 2002, hereafter WPM02). 

Of the 150 NLSls hstcd in WPM02, 52 were detected in the ROSAT All-Sky Survey 
(RASS; Voges et al. 1999). Forty-five of these had sufficient counts in the 0.1-2.0 keV range 
to derive power-law photon indices (F, where N{E) oc E^^) based on hardness ratios. Most 
of these objects were optically bright {g < 18.5), low-redshift {z < 0.4), and exhibited 
typical NLSl photon indices of F > 2.0. However, a substantial number of optically-bright, 
low-redshift NLSls did not have RASS source identifications. The optical spectra of these 
objects appear completely normal (within the limitations of the SDSS resolution and signal- 
to-noise) and only one or two are in regions of high Galactic H I column density, which could 
potentially obscure the X-ray fiux. 

It is thus possible that these objects represent a subset of NLSls which are optically 
normal but X-ray weak. Only a few such objects have previously been found; for example, 
RX J2217.9-5941 (Grupe, Thomas, & Leighly 2001) and PHL 1811 (Leighly et al. 2001). 
Since these SDSS NLSls were not detected in the RASS, however, the nature and extent of 
their X-ray weakness could not be determined. They could emit significant fiux at higher 
energies not covered by the 0.1-2.4 keV ROSAT band, or they might also be ultrasoft X- 
ray sources but with substantially lower overall X-ray fluxes than typically seen in NLSls 
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(i.e., much higher ctox)- Another possible explanation is variability, but it seems unlikely 
that such a large fraction (~ 40%) would be in an exceptionally low state during the RASS 
observations. In reality, all three of these factors probably have some bearing on the X-ray 
weakness of these objects, but we cannot infer how many are affected by which factors, if 
any, from the existing data. 

As a first step toward solving this puzzle, we have observed 17 of these optically- 
selected but RASS-undetected NLSls with Chandra. Due to its excellent sensitivity and low 
background levels, Chandra is able to detect objects at far lower fiux levels than ROSAT, 
and its large (0.5-8 kcV) energy range allows detection of objects with harder X-ray spectra 
as well. Our primary goals are (1) to detect these objects in X-rays, or set upper limits to 
their X-ray emission, and (2) to obtain rough estimates of F for the C/ianc?ra-detected NLSls. 
Given this information, we can gain some insight as to which of the aforementioned scenarios 
(if any) sufficiently explain these ROSAT-unobserved NLSls. In this paper we present the 
results of these Chandra observations, and the possible implications for the NLSls in our 
sample. 



2. Target Selection 

A full description of the spectroscopic selection and preliminary X-ray analysis can 
be found in WPM02; a brief summary follows. The NLSls in the WPM02 sample were 
initially selected from the SDSS Early Data Release (Stoughton et al. 2002) solely on the 
basis of narrow H/3 emission. Subsequently, each spectrum was visually inspected and a 
more accurate measurement was taken of the H/? line width. Objects which did not fulfill 
the criteria of Ostcrbrock & Poggc (1985) and Goodrich (1989) were discarded. X-ray power 
law slopes and luminosities were estimated for the 45 objects with adequate RASS data 
available. 

Figure 1 shows the SDSS g magnitude-redshift distribution for the WPM02 sample, 
differentiating between NLSls detected in RASS, those not detected, and those chosen for 
Chandra follow-up. In this diagram, the fainter, higher-redshift NLSls are detected in the 
RASS less frequently, and most NLSls with g < 18.5 and z < 0.4 were detected by RASS. We 
adopted these as rough limits for our follow-up sample. To rule out the possibility that high 
foreground extinction prevented some objects from being detected in the RASS, we restricted 
our sample to objects in regions of low Galactic H I column density^ (A^h < 4 x 10^° cm~^. 



^ Found for each location using the nh utility, part of the HEAsoft package, available at 
http: / /heasarc. gsfc.nasa.gov/lheasoft/. 



-4- 



although one relatively bright NLSl in the sample has A^h = 5.7 x 10^° cm The resulting 
17 objects which comprise this follow-up sample appear normal in all other respects. 

3. Observations and Data Reduction 

We observed the 17 NLSls in this sample with Chandra for ~ 2 ksec each between 3 
January and 15 September 2003. The Advanced CCD Imaging Spectrometer Spectroscopic 
array (ACIS-S; Garmire 2003) was employed, and all observations were offset by 3' from the 
nominal aimpoint along each detector axis to mitigate the effects of pilcup should any of 
the sources be unexpectedly bright. This resulted in a slight (~ 10%) reduction in detector 
efficiency as well as a broadening of the point-spread function. Neither of these effects 
significantly hinders our source detection efficiency, thanks to Chandra's exceptionally low 
background levels. These observing parameters allow 3a detections of point sources down 
to 2 - 5 X lO-^'^ergcm-^s-^ (for T = 3.0 - 1.0 respectively) in the 0.5-8 keV band. 

3.1. Data prepeiration 

All data were reprocessed with the newest version of the Chandra Interactive Analysis 
of Observations (CIAO 3.0.1) software employing the most recent calibration files (CALDB 
version 2.23).^ We then defined circular extraction regions that were centered on each source 
and of sufficient radius to encompass the PSFs (typically about 4'.'5). Background regions 
were defined as annuli with inner and outer radii of 7'.'5 and 15" respectively. The CIAO 
tool psextract was used to extract source and background spectra for those objects with 
sufficient flux (> 25 total counts). Additionally, we used dmextract to determine the raw 
number of counts for each source in the 0.5-0.9 (S), 0.9-2.0 (M), and 2.0-8.0 keV (H) bands. 
Finally, the apply_acisabs script was applied to each ancillary response file (ARF) in order 
to account for the time-dependent ACIS quantum efficiency degradation.^ 

All 17 NLSls in this sample were detected by Chandra, 16 with high (> Aa) significance. 
Table 1 lists the net, background subtracted, number of counts per energy bin for each object. 
As expected, the background was found to be low, with 0-2 counts per band detected in 
the background extraction regions (i.e. less than 0.25 counts expected in the source region 
on average), for most of the observations. In only three cases were one or more background 



^Both available at littp://cxc. harvard.edu/ciao/downloadjciao_reg. html 
•^See http://cxc.harvard.edu/cal/Acis/CaLprods/qeDeg/ for details. 
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counts expected in the source region. Due to the low background levels, we assume the errors 
are simply Poisson errors on the raw number of source counts. 

3.2. Determination of Power-Law Slope 

We used the CIAO spectral fitting tool Sherpa to analyze the resulting data. Back- 
ground levels are negligible (< 1 count expected in the source extraction region) in all but 
two observations; for these two, we fit the background simultaneously with a power law 
model. For the twelve objects with sufficient counts to obtain binned spectra (with at least 
5 counts per bin), we fit a simple power law with foreground Galactic absorption taken into 
account"^. This model can be represented by the equation 

dN{E) cx E-^e-^"''^^^A{E)dE (1) 

where (t{E) is the photoelectric absorption cross section of Galactic gas with effective column 
density Nh, and A{E) is the observation-specific effective detector area encoded within each 
ancillary response function (ARF) file used within Sherpa. Knowing Nhj T and the overall 
normalization can be found through a simple two-parameter fit to the data. 

For those objects too faint to be analyzed as binned spectra in Sherpa (and to check the 
consistency of our spectral fits) we employed a hardness ratio fitting method. This method 
reduces the problem to a one-parameter fit by disregarding the overall normalization of the 
spectrum (which can be determined afterwards from the total count rate, {H + S)/texp)- 
The hardness ratio is defined as HR = {H — S)/{H + S), where H and S are now the 
net counts in some arbitrary hard and soft bands respectively. Assuming Poisson errors, 
chr = 2\/WlT^\^lPS / {H -\- Sy (again using the raw; number of counts per band to calculate 
errors, for those observations with non-zero background levels). We employed two different 
hardness ratios in this analysis, which we call HRa and HR;,. HRa is analogous to the 
ROSAT "hardness ratio 2" parameter, where the soft and hard bands are 0.5-0.9 and 0.9- 
2 keV respectively. HR5 takes advantage of the full Chandra energy range with soft and hard 
bands covering 0.5-2 and 2-8 keV, respectively. Again using the ARF and Nh value specific 
to each observation, we used Sherpa to calculate HRa and HR;, for a test value of F. By 
iterating F until the correct values of the hardness ratios were reached, we derived photon 
indices independently for each of HR^ and HR^. 



^For fitting we used the statistic with the Gehrels (1986) variance function, which is the Sherpa default 
and applicable to data with few counts. 
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3.3. Consistency 

For 15 of these NLSls, F was derived by using two or three of the methods discussed in 
the previous section. As a consistency check, we compared the resuhs for individual objects. 
For the most part, especially for the brighter objects, the three measurements produced 
results well within la of each other. However, there are systematic offsets. Values of F 
derived from HRa are nearly always higher than those found using HR;,, while Sherpa fits 
typically lie somewhere in between. This is as expected since the soft X-ray power-law slopes 
of NLSls are known to be steeper than those in the hard band due to the commonly-observed 
soft excess (e.g. ROSAT: BBF96 vs. ASCA: Brandt, Mathur & Elvis 1997, although some 
of this offset may be due to calibration issues as described in Iwasawa, Fabian, & Nandra 
1999). A slope over the entire energy range, as measured by Sherpa, would be in between 
the two. Since they contain the average spectral information, we list the F from Sherpa 
fits in Table 2. If there are not enough counts for such a fit, we hst the 0.5-8 keV (HRfe) 
measurement since it covers the full energy range; for the lone observation with no 2-8 keV 
counts, the 0.5-2 keV measurement of F is listed. 

4. X-ray Spectral Properties 

4.1. Detectability in the RASS 

The NLSls in this sample span a broad range of F (best fit values of 0.25-3.15) and 
X-ray fluxes from near the detection limit of about 0.002 counts to almost 0.2 counts s~^. 
Taking the 0.5-2 keV flux derived from Sherpa for the twelve brightest objects, we used the 
PIMMS software^ to calculate the count rate expected in ROSAT. Exposure times for each 
position were taken from the RASS exposure maps^ and the expected number of RASS counts 
was then computed for each object. Three of the sources, J1013-1-0102, J1214+0055, and 
J1449-I-0022 should have been easily detected with about 17, 17, and 22 counts, respectively. 
The second object may actually correspond to a nearby RASS source which was marginally 
detected and not cross-referenced in the SDSS database. Four other objects in this sample 
should be marginally detected by the RASS with ~ 10 counts; their nondetection may be 
due to Poisson error or uncertainty in the flux determination. The remaining ten objects 



^Portable, Interactive Multi-Mission Simulator Version 3.4, from NASA's High 

Energy Astrophysics Science Archive Research Center, currently available at 
http://heasarc.gsfc.nasa.gov/docs/software/tools/pimms.html. 

^Available at http: / /www. xray.mpe.mpg.de/cgi-bin/ rosat / rosat-survey. 
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were well below the RASS detection limit. 

The X-ray luminosity of the three bright objects may have varied by a factor of two or 
more between their RASS and Chandra observations, which could account for their lack of 
detection in the RASS. Such long-timescale variability is not surprising, as it is commonly 
seen in NLSls (e.g., Grupe, Thomas, & Beuermann 2001). It is impossible to ascertain 
the degree of variability in the fainter objects; however, most of this sample could not have 
been detected in the RASS at the flux levels observed by Chandra. It is unlikely that all 
14 of these 17 NLSls were in an exceptionally low luminosity state during both the RASS 
and Chandra observations; thus, at least some of these objects must be intrinsically faint at 
ROSAT energies. 

4.2. Photon Index— Luminosity Relation 

It is likely that some combination of variability, X-ray hardness, and low X-ray lumi- 
nosity resulted in these objects not being detected in the RASS. The latter two factors can 
fortunately be determined from these observations. The data in Table 2 show that many 
of the NLSls with low F are also the faintest X-ray sources. To determine whether this is 
a true effect, we calculated the monochromatic 1 keV (rest- frame) fluxes using the models 
derived through spectral fitting in Sherpa. For the fainter objects we calculated the expected 
flux from the best-flt hardness ratio models using Sherpa, with a fractional flux error equal 
to N'^f^ where is the total number of counts in the 0.5-8 keV band. Using the SDSS 
redshifts and an Vt^ — 0.3, Oa = 0.7, Hq = 70kms~^Mpc~^ cosmology, the corresponding 
rest-frame luminosities were then calculated. These are given in Table 2, along with the 
optical A5100A monochromatic flux measured from the original SDSS spectra. 

The relation between F and Likev is shown in Figure 2. Luminosities derived from 
Sherpa fitting are plotted as filled circles, and those estimated from the hardness ratio 
models are plotted as open circles. This figure shows a strong correlation between X-ray 
spectral slope and 1 keV luminosity. Grupe et al. (1999) found such a correlation among 
soft X-ray selected AGN, with a stronger correlation among NLSls noted in Grupe (2004). 
No such correlation is seen in their broad-line AGN sample; additionally, Laor et al. (1997) 
do not find a signiflcant correlation in their sample of bright, optically-selected quasars 
(although this may be due to the hmited range of F in these samples). 

This relation has a slope of 6 = 0.66 ± 0.13 (where F oc 61ogLikev); if the two hardest 
objects are excluded from the flt, then b decreases to 0.55 ± 0.13. Similar relations have 
been measured for NLSls by Forster & Halpern (1996) with b — 0.32 (where the luminosity 
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is measured over the 0.1-2 keV range), while the Dai et al. (2004) quasar sample has b — 
0.55 ± 0.11 (with L = Lo.2-2kev)- When the luminosities of our sample are extrapolated to 
the 0.1-2 keV range using the best-fit models and the two lowest-F points are excluded, the 
slope becomes b = 0.48 ±0.08, in agreement with Dai et al. (2004) but about 2o" higher than 
the Forster & Halpern (1996) result (note that the choice of 0.1 keV or 0.2 keV as the lower 
bound of the energy range does not significantly affect the extrapolated luminosity). This 
latter discrepancy may be due to their sample being selected from the RASS, which is not 
as sensitive to fiat-spectrum objects as Chandra. 

For contrast, the F — Likev relation caused by pure absorption (i.e. assuming a typical 
source with F ~ 2.5, determining what hardness ratio and Likev would be observed with 
various intrinsic Nh values, and recalculating F using the simulated HR^ assuming no prior 
knowledge of the intrinsic absorption) produces a slope closer to 6 = 2 for Nh < 5 x 10^^ cm~^. 
Most of the observed correlation thus cannot be explained simply by intrinsic absorption, 
although it could explain the hardest objects in this sample (F < 1). However, this would 
require large neutral column densities (> lO^^cm"^), despite there being no sign of optical 
dust reddening. Moreover, it is unlikely that all four NLSls with F < 1 arc strongly absorbed 
(see §4.4). A dust-poor absorber could, in principle, account for the extreme hardness and 
low X-ray luminosity of one or two of these objects. 

The observed F — Likcv correlation may be a result of differing Eddington ratios. To 
test this, we estimated black hole masses and Eddington luminosities using the Kaspi et 
al. (2000) relation between i?BLR, ALa(5100A), and FWHM(H/?) and applying the virial 
theorem. It should be noted that these estimates may be subject to systematic error since 
the Kaspi et al. (2000) sample contains very few NLSls; however, this method should be 
sufficient to show general trends in the data. Values of FWHM(H/3) were taken directly 
from WPM02, and the 5100A flux was measured from the SDSS spectra used in that study. 
Figure 3 demonstrates a strong correlation between F and -C/ikevZ-^'Edd; which further indicates 
that the objects with the softest X-ray spectra have the highest relative accretion rates if 
Likev oc Lbol- This is probably a reasonable assumption. Figure 4 shows that the 1 keV 
luminosity is correlated with the 5100A luminosity for the bright and faint objects in this 
sample (with the notable exception of one point, which has F ~ 0.25 and may be strongly 
absorbed) . We can thus assume that both the optical and X-ray luminosities provide some 
indication of the bolometric luminosity. As expected, a correlation is also seen between F and 
ALa(5100A), though it is somewhat weaker than that between F and -LikeV- This weakness 
is not surprising, since X-ray properties should physically be more closely linked with each 
other than with optical measurements. 

Although F is directly proportional to i^boiZ-^Edd, if is thought to be inversely propor- 
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tional to Mbh (e.g., Kuraszkiewicz et al. 2000), which could wash out the observed F — Likev 
correlation for some samples. For this Chandra sample, the Mbh estimates span approxi- 
mately two orders of magnitude, while Likcv / -^Edd spans three, so the correlation is observed. 
In samples with a larger range of black hole masses and luminosities (such as the Grupe 2004 
BLSls), this correlation would not be expected and, indeed, is not seen. The F — Likcv re- 
lation would also not be seen in samples with a smaller range of ^boi/-^Edd, unless the black 
hole mass range is correspondingly larger, as proposed for the Dai et al. (2004) sample. 
When the same analysis is performed for the WPM02 RASS objects, no correlation is ob- 
served. Instead, they cluster almost exclusively around the highest-Likev/-^Edd, highest-F 
members of the Chandra sample. This is not unexpected since the WPM02 RASS sample 
consists mostly of NLSls near the RASS detection limit; they do not exhibit the range in 
luminosities (and hence i^ikcv/-^Edd) or required to sec the F — Likcv correlation. If the 
observed F — Likcv correlation is indeed due to differences in the Eddington ratio, then it 
is quite likely that this Chandra sample includes some NLSls which are accreting far below 
the Eddington limit. 

4.3. Spectral Features in J1449-|-0022 

As noted in §4.1, J1449+0022 is the brightest object (in X-rays) in this sample and may 
have gone undetected in the RASS due to strong variability. Indeed, there appears to be a 
slight decrease in the Chandra count rate over the course of the 2 ksec observation. More 
interesting, however, is that the spectrum of this object is not well fit by a simple power law 
with Galactic foreground absorption, as seen in Figure 5; specifically, there appears to be 
an excess of soft photons over that expected from a power law. The X-ray continuum can 
be better fit using a slightly more complicated model; either a power law plus a blackbody 
(which yields F = 1.58 ± 0.22 and kT = 0.11 ± 0.02), or a broken power law with Fi = 
SAOtoll, F2 = 1.64 ± 0.19, and Eb = 1.041°;}^ kcV. There is also a significant dip in fiux at 
1.10±0.03keV with equivalent width 0.121qq4 keV. This absorption is similar to that seen in 
some other NLSls (Leighly et al. 1997). Such a feature could be due OVII/OVIII absorption 
in a highly relativistic outfiow (e.g.. Pounds et al. 2003), but a more plausible model may 
be the Nicastro, Fiore, & Matt (1999) hypothesis of a strong Fe L complex at this energy. 

4.4. Intrinsic Absorption 

Four of the NLSls in this sample exhibit spectra with F < 2. Since for low count rates 
F and A''^^ are highly degenerate, no individual object is bright enough to place a meaningful 
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limit on its intrinsic absorption; thus, we cannot determine with confidence whether this 
X-ray hardness is intrinsic to the accretion process or merely caused by a high degree of 
absorption. The spectra of these four NLSls were added together to determine whether 
they exhibit strong absorption as a group. We assumed all have typical T = 2.5 and similar 
Nh values at each redshift, and fit a corresponding model to the coaddcd spectrum. This 
resulted in a best-fit value of Nh = l-Slgig ^ lO^^cm"^, but the fit is not particularly good 
and shows strong residuals. If we leave both F and Nh as free parameters, Nh becomes zero 
(with a 2(7 upper limit of 2 x lO^^cm"^) and F = 1.2 ±0.3. The fit with no absorption is also 
much better (Ax^ = —6.6 compared to the absorption model). Based on this, it appears as 
though these flat-spectrum sources are not strongly absorbed as a group. 

Of course, it is possible that one or two of these NLSls could be strongly absorbed, 
particularly the faintest object (J1259-I-0102) since with only flve total counts its influence 
on the coadded spectrum is small. If this object has an intrinsic X-ray luminosity Likev ~ 
lO^^ergs"^ and F ~ 3 (similar to the brightest objects in our sample), an intrinsic neutral 
column density of ~ 2.5 x lO^^cm"^ would be required to reproduce the observed flux and 
F. The input F and luminosity of such a model can be adjusted to reproduce the other 
faint, hard sources as well, but somewhat high column densities {Nh ^ 5 x lO^^cm^^) are 
typically required. Moreover, there is no indication of absorption in the optical spectra of 
these objects. Thus, we conclude that these objects most hkely have intrinsically hard X-ray 
spectra, although the two hardest may be heavily obscured in the X-rays but not at optical 
wavelengths (cf. Risaliti et al. 2001; Nandra et al. 2004). 



5. Compcirison to Other NLSl Samples 

5.1. X-ray Spectral Slope 

At first glance, this sample appears to exhibit a range of X-ray spectral slopes extending 
to much lower values than previously seen for NLSls. However, the F values listed are, for the 
most part, derived over a larger energy range than used by previous studies. For example, F 
is measured by ROSAT over the 0.1-2.4 kcV energy range for the 52 narrow- and broad-line 
Seyfert 1 objects in BBF96, while the 0.1-2.0 keV hardness ratio is used with the WPM02 
sample. In order to ensure the energy bands overlap, we employ the HRq measurement of F 
for comparison to these two samples. 

It should be noted that the energy ranges used by these samples arc slightly different 
nonetheless. Since HR^ covers the 0.5-2.0 keV band, it docs not extend to the 0.1 kcV 
minimum energy of ROSAT, or to the 2.4 keV maximum energy used by BBF96. However, 
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28 of the WPM02 NLSls detected in ROSAT have measurements for both HRl (0.1-2.0 keV) 
and HR2 (0.5-2.0 keV). We rederived F using the HR2 measurements and compared it to the 
values from HRl. Although the errors are larger, the HR2 Fs are consistent with those from 
HRl in all but four objects. No systematic offset is seen between these two measurements. 
Thus, we assume that the Fs derived from ROSAT HRl are good estimators of the spectral 
slopes over the 0.5-2 keV range. As for the BBF96 measurements extending to 2.4 keV 
while ours extend only to 2.0 keV, we expect this to be neghgible for two reasons: first, few 
photons are emitted in the high-energy tail (A^2.o-2.4keV ~ 0.1A^o.5-o.9keV for T = 2); and 
second, the ROSAT sensitivity decreases quite rapidly in this regime, so even fewer photons 
should be detected. Thus, the differing energy bands should have httle effect on comparing 
these various samples. 

Figure 6 shows the relation between F and the H/3 velocity width for all three samples. 
Crosses denote BBF96 data, blue circles are the WPM02 RASS data, and the 16 Chandra 
objects with HR^ measurements are shown as red circles with errorbars. The horizontal lines 
show the mean and sample standard deviation of F for only the BBF96 data; NLSl and Syl 
averages are computed separately. It is immediately apparent that the WPM02 and Chandra 
NLSls do not fall within the same range as those in BBF96; in fact, a significant number 
appear to exhibit Seyfert i-type X-ray spectra. Table 3 lists the mean and sample standard 
deviation of F for the three samples. An interesting progression emerges: as we move from 
NLSls which were found through a mixture of optical and X-ray selection (BBF96), to 
those which were selected solely based on their optical properties ( WPM02) , to the Chandra 
sample presented herein which was selected on the basis of weak X-ray emission, (F) and 
the extremes of the F distribution become harder. In fact, the F distribution of this sample 
extends to lower values than that of the BBF96 Seyfert 1 galaxies. This illustrates the effect 
of selection methods on building NLSl samples, as well as demonstrating the existence of a 
significant population of NLSls which are not ultrasoft X-ray sources. 



5.2. X-ray vs. Optical Luminosity 

This Chandra sample was originally selected as a set of X-ray weak NLSl candidates, 
i.e. NLSls which were not detected in the RASS but had g' magnitudes similar to detected 
objects. To determine whether these objects are truly X-ray weak, we used the WPM02 
RASS-detected NLSls as a comparison sample. Rest-frame luminosities at IkeV were 
calculated using PIMMS along with ROSAT count rates and the best-fit F values given in 
WPM02, and 5100 A luminosities were measured from the original SDSS spectra. The best- 
fit relation between the monochromatic optical and X-ray luminosities has a slope consistent 
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with unity, and is plotted on Figure 4 along with the data points from the Chandra sample. 

Interestingly, the RASS relation corresponds very well with the upper bound on the 
X-ray-optical relation for the Chandra NLSls. Eight of these objects are consistent with the 
RASS fit and another three are consistent within ~ 1 — 2(7, while the rest are significantly 
below the line. The faintest objects in X-rays appear to exhibit a Lx — i^opt slope similar 
to that of the brighter objects (if the point at the lower right of the plot, which is probably 
absorbed in X-rays, is excluded), but their X-ray luminosities are roughly a factor of 5 lower 
than those seen in the WPM02 RASS sample. Thus, although most of the objects in this 
sample are not particularly X-ray faint compared to NLSls in the WPM02 RASS sample 
with similar optical luminosities, six of them do fall well below the Lx — i>opt relation due 
to intrinsic X-ray faintness or obscuration. 



6. Discussion and Conclusions 

Through short-duration Chandra observations of RASS-undetected NLSls, we have 
determined that six of the 17 objects have X-ray luminosities substantially lower than NLSls 
with similar optical properties. Of the brighter objects, at least two exhibit flux levels 
which should have been detectable by the RASS, indicating that their luminosities may have 
increased by a factor of two or more between the RASS and Chandra observations. Many 
of the remaining bright objects were near or just below the RASS detection limit, and were 
most likely not seen due to smaller luminosity variations or Poisson noise. Across the entire 
sample, a strong correlation is seen between the X-ray spectral slope F and Likev- This is 
probably not entirely due to intrinsic absorption, since individual spectra of bright objects as 
well as a coadded spectrum of the faintest objects do not indicate high degrees of absorption 
(although one or two of the faintest hard-spectrum objects may be absorbed in X-rays but 
not at optical wavelengths). If F is indeed correlated with Lboi/-^^Edd, then the F — Likev 
relation suggests that variations in Likev are primarly due to differences in Lboi/ -^Edd among 
objects with comparatively similar black hole masses. This interpretation is complementary 
to that of Dai et al. (2004), who find a similar relation but whose sample more likely includes 
objects with a large range of Mbh but I/boi/-t'Edd ~ 1- 

These observations may hold important implications for the "Eigenvector 1" (Principal 
Component 1; PCI) paradigm posited by Boroson & Green (1992) and reinforced by Boroson 
(2002). In this picture, PCI (which is primarily driven by an anticorrelation between [O III] 
and Fe II) is an indicator of Lboi/-Z^Edd- NLSls typically lie at one extreme end of PCI — the 
end thought to correspond to the highest relative accretion rates. Since F is thought to be 
related to ^boi/-^Edd) T and PCI should be correlated; indeed, Brandt & Boiler (1998) find 
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such a correlation. However, the sample presented herein contains objects which from their 
optical spectra are at the supposed high-Lboi/-^Edd end of PCI, yet also exhibit very low 
values of F, as well as low inferred -^ikev/-^Edd- 

These extreme objects may indicate that while PCI is usually correlated with I/boi/-^Edd, 
it may also be affected by orientation, black hole mass, or other physical drivers (as noted 
by Boroson 2004). This is not a completely new phenomenon; for example, BBF96 note 
that Mrk 507, with F = 1.6 ± 0.3, has an unusually flat X-ray spectrum for a NLSl. This 
sample simply demonstrates that Mrk 507 is not an isolated case, and in fact a small but 
interesting subset of NLSls do not appear to fit within the PCI framework. The apparent 
lack of strong absorption in some of these flat sources indicates that their X-ray spectra 
actually are intrinsically flat. Further studies of X-ray weak NLSls, as well as much larger 
samples from surveys such as the SDSS, should offer greater insight into the mechanism(s) 
behind PCI. 

Due to the short exposure times of these observations (< 2 ksec), we cannot infer 
much outside of F and luminosity estimates for individual objects; indeed, this program was 
intended to study the group properties of an X-ray weak NLSl sample. However, there are 
several objects in this sample with exceptionally low F which may be worthy of further study. 
These hard X-ray NLSls may represent a new, rare subclass which are optically normal but 
highly absorbed in the X-rays, or which exhibit abnormally low i^boiZ-^Edd, or both. 

The authors are grateful to Th. Boiler for providing the original BBF96 data for our 
Figure 6. We also thank the anonymous referee for helpful comments and suggestions. 
Support for this work was provided by NASA through Chandra Award Number G03-4145X 
issued by the Chandra X-ray Observatory Center, which is operated by SAO for NASA under 
contract NAS8-39073. 
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Fig. 1.— SDSS g' magnitud^redshift relation for the WPM02 SDSS-selected NLSls. 

Crosses indicate objects undetected in the RASS, open circles are RASS-detected objects, 
and filled circles are NLSls undetected in the RASS which were chosen for C/ianc/ra follow-up. 
The dashed lines show the nominal brightness and redshift limits for the Chandra follow-up 
sample. 
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Fig. 2. — r (from Table 2) vs. 1 keV rest-frame luminosity for the Chandra foWow-up NLSl 
sample. Filled points are bright objects with spectra fit in Sherpa, while for open circles the 
luminosity was estimated using HR^ (HRa for J131 1+0031) and the 0.5-8 keV count rate. 
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Fig. 3. — Same as Figure 2, but with photon index plotted against -Likev/-^Edd- The x- 
axis errorbars only reflect the uncertainty on Likev since the uncertainty in LEdd is not 
well-known. 
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Fig. 4. — X-ray vs. optical luminosity for the objects with F derived in Sherpa (solid points), 
or using hardness ratios (open circles). The solid line denotes the best-fit LikeV— -^opt relation 
fi:om the WPM02 RASS sample. The high-optical, low-X-ray luminosity point is the object 
with F = 0.25, and may be highly absorbed. 



-20- 




1 

Energy (keV) 



Fig. 5. — Chandra count spectrum for J1449+0022, along with the best-fit power-law model 
and residuals. This spectrum has been binned to 10 counts per bin. The observed flux is 
systematically high at £^ < O.SkeV, and there is an absorption feature just above 1 keV. 
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(j) xapui uo-;oqd ^bj-x Ijog 

Fig. 6. — Soft X-ray photon index vs. H/? width, for NLSls and Seyfert 1 galaxies (on 
the left and right of the dashed vertical line, respectively). Crosses are data from BBF96, 
kindly provided by Th. Boiler. Blue circles are NLSls from WPM02 (with F derived from 
RASS measurements) and red circles are the sample presented herein. All values of F for 
the Chandra follow-up sample were found using the 0.5-2 keV hardness ratio HRq. The 
horizontal solid and dashed lines show the means and sample standard deviations of F for 
the BBF96 NLSls and Syls. 
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Table 1. Chandra observing log 



SDSS Name^ 




Date 


^exp 

(s) 








(1020 cm-2) 


is-') 


J002305.03-010743.5 


25- 


-08- 


-2003 


1940 


30 


44 


13 


2.81 


0.045 ± 0.005 


J002752.39+002615.8 


03- 


-09- 


-2003 


1560 


7 


21 


14 


2.72 


0.027 ±0.004 


J015652.43-001222.0 


12- 


-09- 


-2003 


1710 


8 


6 


3 


2.62 


0.010 ±0.002 


J022756.28+005733.1 


23- 


-06- 


-2003 


1730 


6 


10 


7 


2.84 


0.013 + 0.003 


JU22841.48+UU52U8.0 


24- 


-06- 


-2003 


OA A A 

2000 


51 


65 


22 


2.7b 


A A/^A 1 A AA/^ 

0.069 + 0.006 


J031427.47-011152.4 


03- 


-09- 


-2003 


1910 


48 


43 


9 


5.70 


0.052 + 0.005 


J101314.86-005233.5 


08- 


-01- 


-2003 


1990 


63 


70 


14 


3.64 


0.074 + 0.006 


J104230. 14+010223.7 


20- 


-02- 


-2003 


1730 


33 


42 


12 


3.72 


0.050 + 0.005 


J121415. 17+005511.4 


07- 


-02- 


-2003 


1960 


49 


42 


12 


1.94 


0.053 ± 0.005 


J125943.59+010255.1 


04- 


-03- 


-2003 


1940 





2 


3 


1.62 


0.003 ± 0.001 


J131108.48+003151.8 


10- 


-03- 


-2003 


1730 


7 


6 





1.90 


0.008 ± 0.002 


J141234.68-003500.0 


07- 


-01- 


-2003 


2090 


36 


48 


16 


3.29 


0.048 ± 0.005 


J143030.22-001115.1 


23- 


-04- 


-2003 


1950 


5 


10 


9 


3.15 


0.012 ±0.003 


J144932. 70+002236.3 


09- 


-07- 


-2003 


2150 


129 


165 


68 


3.75 


0.168 ±0.009 


J145123.02-000625.9 


22- 


-04- 


-2003 


2120 


52 


74 


14 


3.84 


0.066 ±0.006 


J170546.91+631059.1 


17- 


-09- 


-2003 


1940 


6 


8 


2 


2.57 


0.008 ±0.002 


J233853.83+004812.4 


28- 


-08- 


-2003 


1910 


28 


23 


11 


3.88 


0.032 ± 0.004 



^Format: SDSS JHHMMSS.ss±DDMMSS.s 

•^Net counts in soft (0.5-0.9 keV), medium (0.9-2.0 keV), hard (2.0-8.0 keV) bands. 
^Galactic foreground H I column density, found with the nh utility. 
'^Net Chandra 0.5-8 keV count rate. 
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Table 2. NLSl Optical and X-ray Properties 



iName 




z 


r VV riiVH^rip ) 
(kms ) 


J0023-0107 





166 


1160 


J0027+0026 





205 


1830 


J0156-0012 





163 


1320 


J0227+0057 





128 


770 


J0228+0052 





186 


990 


J0314-0111 





387 


1810 




n 


976 


1 580 


J1042+0102 





116 


1010 


J1214+0055 





396 


1980 


J1259+0102 





394 


1460 


J1311+0031 





429 


1640 


J1412-0035 





127 


1100 


J1430-0011 





103 


1740 


J1449+0022 





081 


1070 


J1451-0006 





139 


1120 


J1705+6310 





119 


1660 


J2338+0048 





170 


1010 



\og{ELEY^a \og{\L^Y'nog{^f Note 
(ergs-^) (ergs-^) 



2.21 ±0.41 


42 


89 


± 





20 


43 


89 


6 


73 


1.61 ± 0.50 


42 


80 


± 





29 


43 


90 


7 


14 


2.04+°-I^ 


42 


24 


± 





24 


43 


84 


6 


82 




49 


04 




n 


91 


43 


40 


U 


04 


2.46 ±0.25 


43 


27 


± 





09 


43 


85 


6 


57 


3.15 ±0.44 


43 


97 


± 





16 


44 


70 


7 


69 


2.71 ±0.34 


43 


74 


± 





08 


44 


43 


7 


38 


2.38 ±0.31 


42 


67 


± 





14 


43 


66 


6 


46 


2.70 ±0.28 


43 


93 


± 





13 


44 


50 


7 


62 


0.25t?:^? 


41 


85 


± 





45 


44 


54 


7 


39 


3.06 ±0.96 


43 


27 


± 





28 


44 


70 


7 


60 


2.34 ±0.38 


42 


74 


± 





14 


43 


36 


6 


31 


0.92 ±0.64 


41 


62 


± 





47 


43 


19 


6 


60 


2.12 ±0.16 


42 


82 


± 





08 


43 


38 


6 


31 


2.56 ±0.22 


42 


96 


± 





12 


43 


57 


6 


48 


2 40+^-°2 

^•^U_Q 55 


41 


96 


± 





25 


43 


55 


6 


81 


2.69 ±0.69 


42 


86 


± 





18 


43 


63 


6 


43 



1 

2 



3 
1 



^Truncated to JHHMM±DDMM 

''Redshift from the SDSS catalog, as listed in WPM02 

•^From WPM02 

•^Derived from spectral fitting in Sherpa, unless noted otherwise. 

^Monochromatic, rest-frame 1 keV luminosity inferred from the best-fit Sherpa or HRa/b model 
fiux, with quoted errors from Sherpa or Poisson errors on the total number of counts respectively. 

^Monochromatic, rest-frame 5100A luminosity. Fluxes are measured from the original SDSS 
spectra used in WPM02; errors are considered neghgible for Mbh estimates. 
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^AU luminosities are calculated assuming Hq — 70kms ^ Mpc ^, flm — 0.3, JIa = 0.7. 
''Estimated from the Kaspi et al. (2000) relation. 

Note. — (1) r derived from HRb; (2) T derived from HR^; (3) Spectrum not well fit with power 
law + Galactic absorption model; see §4.3. 
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Table 3. Statistics of T for NLSl X-ray 
Samples 



Sample 


(r) 


CTr 


r Range 


BBF96 NLSls 


3.22 


0.73 


1.8- 


-5.2 


WPM02 (RASS) 


2.75 


0.31 


1.5- 


-4.4 


This Paper 


2.60 


0.39 


1.1- 


-3.4 


BBF96 Syls 


2.31 


0.28 


1.4- 


-3.4 



